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Abstract 

Recently, a novel framework to handle stochastic processes emerges from a series of studies 
in biology, showing situations beyond "Ito vs. Stratonovich" . Its internal consistency can be 



Cd demonstrated via the zero mass limit of a generalized Klein-Kramers equation. Moreover, the 

^ connection to other integrations becomes evident: The obtained Fokker-Planck equation defines a 

H new type of stochastic calculus that in general differs from the a-type interpretation. A unique 

'^ advantage of this new approach is a natural correspondence between stochastic and deterministic 

^ dynamics, which is useful or may even be essential in practice. The core of the framework is 
a transformation from a usual Langevin equation to a form that contains a potential function 

I „. .„ ..... .„ .... ... ™. . ..« ,„ „. .. 

^^ framework has a direct physical meaning and a straightforward experimental realization. A recent 



experiment offers a first empirical validation of such new stochastic integration. [This article has 



O been published in J. Stat. Meek. (2012) P07010[ Note that a typo in Sect. Ill B has been corrected 



^ (See also [39]). ] 
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I. INTRODUCTION 

Ito-Stratonovich dilemma [IH3], on choosing the appropriate calculus when integrating 
stochastic differential equations (SDEs) "attracted considerable attention in the physics 
community" and "is still as elusive as ever" [3] . Recent explorations of SDEs [3HTU] have well 
been expanded into biology and other fields, generating new understandings on stochastic 
calculus. The novel framework [GlIOHTT] in the present brief review is beyond Ito-Stratonovich 
controversy from two aspects: First, it defines a new type (named A-type for short) of 
stochastic integration that is different from Ito and Stratonovich's. In one dimensional case, 
A-type integration reduces to the a-type |1] with a = 1 (Ito's corresponds to a = and 
Stratonovich's to a = 0.5). In higher dimensional situations, it is generally not an a-type 
integration ^E\. Second, there is a starting point in which obtaining this new approach, we 
need not to choose the "correct" calculus. One can reach the same result by Ito, Stratonovich 
or A-type prescription of stochastic integration. 

The effectiveness of the new framework has been demonstrated by its various applica- 
tions: The successful solution of the outstanding stability puzzle of a genetic switch [5j; the 
quantification of evolutionary dynamics of Darwin and Wallace |7j; the study of complex 
bio-networks such as metabolic network [8J and cancer network [12]; the implications of Dar- 
winian dynamics in physics |13j; the relationship with Lyapunov's direct method for stability 
analysis in engineering [T7]; the explicit construction of adaptive landscape in population 
genetics [15] and etc. Results of a recent experiment in one dimension [18] provide evidence 
that there exist processes in Nature choosing A-type integration. 

A transformation from the classical Langevin equation to a structured form lies at the 
core of this framework. The significance of such transformation is the obtaining of a po- 
tential function that plays a dual role: It leads to the Boltzmann-Gibbs distribution on the 
final steady state (if it exists) of the stochastic process; it corresponds to the determin- 
istic dynamics as a global Lyapunov function and can be used for stability analysis [T7] . 
The potential function exists for general non-equilibrium processes without detailed balance 
[HI [ini [U], which are usually difficult to handle by previous methods [19J. The behaviors 
near a fixed point, stable or not, have been exhaustively studied [9]; explicit construction 
for a simple limit cycle dynamics has been provided [20] as well. The structured form itself 
reveals the symplectic structure embedded in stochastic differential equations and has an in- 



variant property [TT] during the transformation. Other apphcations based on the concept of 
potential function have also been studied in biology [T11I2IH2S], physics I2nH2H] and control 
theory ^W\- 

In the following, we first review the framework in Sect. Ill] about the transformation and 



its consistency in mathematics. In Sect. Ill, several related issues have been discussed, 



including a straightforward implication on the physical meaning. In particular, the one 



dimensional experiment serving as an example is analyzed in Sect. IIIC We summarize in 
Sect. lYl 

II. THE FRAMEWORK 

A heuristic demonstration of this framework is in the following subsection. It is first 
proposed by one of the authors in [6] . A more rigorous exposition [10] is then provided in 
Sect. HTBl 

A. Review of the Transformation 

Stochastic differential equation, or the Langevin equation in physics, is usually a more 
precise description of reality than the purely deterministic one [TH3]. Additive noise is 
frequently considered for physical systems, but for biological or other complex systems, a 
state variable dependent (multiplicative) noise is often encountered. In this brief review, we 
use the physicists' notation with a multiplicative noise: 

q = f(q) + iV(q)e(t), (1) 

where q, f are n-dimensional vectors and f is a nonlinear function of the state variable q. 
The noise ^(t) is fc-dimensional Gaussian white with the zero mean, (^(t)) = 0, and the 
covariance (^(t)^^(t')) = 6{t — t')Ik- The superscript r denotes the transpose of a matrix, 
5{t — t') is the Dirac delta function, (...) indicates the average over noise distribution, I^ is 
the fc-dimensional identity matrix. The element of the nx k matrix A^(q) can be a nonlinear 
function of q. This matrix is further described by: 

N{q)N-{q) = 2eD{q) , (2) 



where e is a constant quantifying the noise strength and D{q) is an n x n positive semi- 
definite diffusion matrix. Note that if the noise has less than n independent components, 
k < n, -D(q) has zero eigenvalue (s). 

During the study of a biological switch ^, a structured form equivalent to Eq. ([T| was 
discovered: 

[5(q) + A(q)]q = -V0(q) + iV(q)e(t) , (3) 

where S{q) is symmetric and positive semi-definite, A{q) is antisymmetric, 0(q) is the 
potential function, ^(t) is identical to that in Eq. ([I]) and the matrix iV(q) is constrained by 



(see also Sect. Ill A): 




iV(q)iV-(q) = 2eS(q) . (4) 

We note that when S = 0, Eq. d3| has the same structure as the Hamiltonian equation 
in physics. The (n = 2/)-dimensional system (/ generalized coordinates and / generalized 
momentums combine the n-dimensional system) with 



A 



and (f) = H (the Hamiltonian), shows the embedded symplectic structure. 

There are two assumptions implied in this framework: The first one is the equivalence 
of Eq. ([I| and Eq. (pi), that is, they are describing the same dynamical process. From 
Eq. ([3]) to Eq. ([I| is straightforward, once [5'(q) + v4(q)] is nonsingular, which holds when 
the components qi{i = 1, . . . ,n) oi the state variable q = (gi, . . . , g^)"^ are independent. The 
reverse, however, is much more difficult, including to obtain S'(q), ^(q) and 0(q) for general 
dynamics. Nevertheless, we can still assume, without a rigorous mathematical proof but 
rigorous enough from a working scientist's view, the reverse part is held based on a case by 
case construction, that is, considering it as a protocol not a theorem. Hence, we may replace 
q with the right hand side of Eq. ([I]) in Eq. ^ : 

[S{q) + A{q)] [f (q) + iV(q)e(t)] = -V0(q) + iV(q)e(t) . (5) 

The second assumption is the separated equality of the deterministic and stochastic 



dynamics in Eq. ([S]): 

[5(q) + A(q)] f (q) = -V0(q) , (6) 

[S{q) + A{q)] iV(q) = N{q) . (7) 

Intuitively, this assumption on separation is plausible for two different reasons: First, the 
zero mean noise function is nowhere differentiable but the deterministic forces are usually 
smooth (C°°) functions, hence two very different mathematical objects; second, the stochas- 
tic force and the deterministic forces describe different timescales of a physical system with 
different physical origins. By replacing Eq. ^ with an equivalent form (IS]) and plugging 
Eq. ([7]) and Eq. ^ into Eq. Q, we obtain the potential condition Eq. ([s]) and the generalized 
Einstein relation Eq. (jol): 

Vx{[^(q) + A(q)]f(q)} = 0, (8) 

[5(q) + A(q)]D(q)[^(q) - A(q)] = 5(q) . (9) 

In principle, the potential function 0(q) can be derived analytically by solving the n{n — l)/2 
partial differential equations (under proper boundary conditions) Eq. ^, together with the 
n(n + l)/2 equations given by Eq. (|9| (ra^ unknowns in [^(q) + ^(q)] and n^ equations). It 
can also be calculated numerically through a gradient expansion [6]. Detailed results near 
fixed points can be found in ^. Explicit constructions of potential functions for limit cycles 
are contained in [TTl [20] . 

In one dimensional case, A = 0, let e = ksT, if the friction 7 is a constant, then 
S = 'j/kBT, Eq. ^ reduces to SD = 'yD /ksT = 1, namely, the product of the friction 
and diffusion coefficients is a constant, discovered by Einstein [31] one hundred years ago. 
Equation ([o]) is a generalized form of the Einstein relation in two ways: The diffusion matrix 
can be nonlinear dependent of the state variable (a verification experiment has been reported 



) and the detailed balance condition can be broken (v4(q) 7^ 0, see also Sect. II B 2). 
From the generalized Einstein relation Eq. (|9]), we know the symmetric part of [^(q) + 
A{q)]-\ D(q) = 1/2 [{S{q) + A{q))-^ + {{S{q) + A{q))-^y], is the diffusion matrix defined 
in Eq. ^. Hence we have [S{q) + A{q)][D{q) + Qiq)] = / (/ is the identity matrix, 
singularity of -D(q), >S'(q) or A{q) can usually be tolerated) where Q{q) is antisymmetric. 
Rewrite Eq. (|6| as 

f(q) = -[D(q) + Q(q)]V</.(q), (10) 



and by multiplying V0(q) on both sides of Eq. ( 10 ), we reach the Hamilton- Jacobi equation: 



f (q) ■ V</)(q) + V</)(q) ■ D(q) ■ V</)(q) = . (11) 

Finally, by the knowledge from statistical physics, we can consider the potential function 
leads to a Boltzmann-Gibbs distribution 

Mq)=2i;yexp(-^). (12) 

as the steady state distribution (if it exists) for the stochastic process described by Eq. ([3]): e 
as the "temperature" and 0(q) as the "Hamiltonian" . Experimental verification is available 
for some cases, e.g. [18] . 

Hence, the two assumptions to implement the equivalence between Eqs. ([I]) and (|3| lead 
to a potential function for general Langevin equation and to a corresponding Boltzmann- 
Gibbs distribution for the final steady state independent of both matrices A{q) and >S'(q). 
While such distribution function seems special, it should be pointed out that it is one of the 
most successful equations in physics, even broader than that of the stationary Schrodinger 
equation. The framework can be applied for systems with an arbitrary noise strength: 
Eqs. (Is]) and (|9| do not contain e. Particularly, deterministic systems can be regarded as the 
weak noise limit (e — > 0) of stochastic systems. The framework is also able to be extended 
directly to include parameter dependent potential and corresponding dynamical components 
[131 EH], e.g. time t and temperature e dependent 0(q, t,e). We note that the Eqs. ^ and 
(|9| have a close connection to Lyapunov's direct method in engineering [17] . These evidences 
make the assumptions more plausible and appealing. 

This construction further leads to a stochastic integration (A-type) different from con- 
ventional ones such as Ito or Stratonovich. One would ask the mathematical consistency in 
this new approach. We point out that the consistency of A-type integration may be seen 
from two aspects: First, Eq. ([I]) or ([s]) by themselves are not complete descriptions of a 
stochastic process. A proper stochastic interpretation is needed, hence A-type is possible in 
principle. Second, we will show in the next subsection that there is an important situation 
in which one starts without the need to differentiate various known stochastic integrations, 
and reaches the consequence that A-type interpretation is the natural choice for the SDE 
Eqs. (§ and Q. 



B. Demonstration of Consistency 

In this subsection, we summarize the main steps to demonstrate that A-type integration is 
mathematically consistent. The rationality is as follows: The starting point is the important 
stochastic differential equation in physics, Newton equation in 2n dimensions with noise 
and friction. Its corresponding 2n-dimensional Fokker-Planck equation (FPE) is uniquely 
determined, known as generalized Klein-Kramers equation. The valuable feature here is, 
all approaches to stochastic integration discussed in the present brief review lead to the 
same 2n- dimensional FPE. When taking the zero mass limit, the 2n dimensional Newton 
equation reduces to an n dimensional equation exactly in the same form as Eq. ([3]), while the 
corresponding n dimensional FPE derived from the generalized Klein-Kramers equation has 
the ability to permit the Boltzmann-Gibbs distribution as its stationary solution. After this 
zero-mass limit, the connection between the resulting n dimensional FPE and the stochastic 
differential equation is the A-type integration. 

1. Generalized Klein-Kramers Equation 

We begin with a general situation via doubling degrees of freedom, the 2n-dimensional 
equation with noise, proposed in [6]: 

q=- 

13) 

p = -[5(q) + A(q)]^ - Vq0(q) + iV(q)e(t) 

where S'(q), A(q), 0(q), N{q) and ^(t) are identical to those in Eq. ([s]), the subscript of 



Vq means it operates on q only, the same as in Eq. rt3|. We note that Eq. (13) has two 



properties: Firsts by taking the zero mass limit ?7i — )■ 0, p — )■ 0, we can recover Eq. ^ 



from Eq. (13); Second, we wish to emphasize that starting with Eq. (13), we can find that 
different stochastic interpretations Ito, Stratonovich, or A-type, will reach the same result, 
since Vp ■ [S(q) + A(q)] = 0. 

The probability density function in the (q, p) phase space can be calculated through path 
integral as: 

p(q, p, t) = (5 (q - q (t, {^(t)})) 5 (p - p (t, {^(t)}))) , (14) 

where q(t, {^(t)}) and p(t, {^(t)}) are the solution of the ordinary differential equation 



(ODE) reduced from Eq. (13) under a fixed sample path {^(t)} of the sample space of ^{t) 



(a given noise configuration). By taking partial derivative of time on botli sides of Eq. (14), 
we obtain 

5tP = ^Vq5 (q- q) ■ ^5 (p-p)+ ,5 (q-q)Vp5(p-p)-^^ 

= _ I <5 (q _ q) 5 (p - p) { Vq ■ [4^^(10)] + Vp ■ [pPd^)] } 40 
= -Vq ■ ^p - Vp ■ {- [^(q) + A(q)] • P p - Vq0(q)p 

+iV(q)(e(t)(5(q-q),5(p-p))| , (15) 



in whicli p, q and p have the same meaning as in Eq. (14), P{{^}) is the probabihty 



density function for the given noise configuration {^{t)}. From step 1 to step 2 of Eq. (15), 
integration by parts is used. According to the Novikov identity jSBj, for a given functional 
/[e(t),t] of noise e(t). 



_/ smt),t] 
\ mt) 



dt' 



and using 



f^at')dt' 



mt) 



noting that the solution of Eq. (13) can be expressed as 



q(t) = q(0) + / P rft' 
Jo m 



P 



p(t) = p(0) + j^\- [5(q) + A{q)] ^ - Vq0(q) + iV(q)e(t') \ dt' 



then we obtain the relations: 



^q(t) _ 

sat) 2 



iV-(q) 



The last left term in Eq. (15) can be handled as 



(16) 



(17) 



(18) 
(19) 



(20) 

(21) 



A^(q) (mS (q - q) 5 (p - p)) = -2eS{q) ■ -VpP 



(22) 



We finally reach the generalized version of the Klein-Kramers equation [2] : 

dtp = Vp ■ {[A(q) + S{q)] ■ ^ + Vq0(q) + e5(q) . Vpj p - ^ ■ V^p . (23) 



Note that Eq. (23) has a solution (if the partition function Z{e) is finite) 



Mq,p) = ^exp[-^^^±^| . (24) 



Equation (24) is not explicitly time dependent, dtPs = 0, it is a steady state distribution. 
The state variables q and p are separated, which simplifies the process of taking zero mass 
limit to be discussed below. 

2. Zero Mass Limit 



We have mentioned that after taking zero mass limit tti — )■ 0, Eq. (13) reduces to Eq. (pj). 
Therefore, by acting the same limit process on Eq. (23), the corresponding FPE to Eq. ^ 
can be obtained [TO] : 

dtp{q, t) = V ■ [D{q) + Q(q)] ■ [eV + V</)(q)]p(q, t) , (25) 



where D(q) is the diffusion matrix, Q(q) is defined in Eq. ( 10 ), V = Vq, and the derivative V 
in V0(q) does not operate on p(q, t). This FPE has a drift velocity — [D{q) + Q{q)] V0(q) = 
f (q) and a generalized symmetric form V ■ [-D(q) + Q(q)] ■ eVp(q, t). The commonly defined 
probability current density j (q, t) = (ji(q,t), . . . , jn(q,^))^ is: 

J.(q, t) = (/. + A/,)(q)p(q, t) - d, [eA,(q)p(q, t)] (26) 



where f (q) = (/i(q), . . . , /n(q))^, A/i(q) = edj[Dij{q) + Qij{q)] (see also Sect. |IIIB| ), with 
i = 1, . . . ,n and j = 1, . . . , n, -Dij(q) the ith-row and jth-column element of the diffusion 
matrix -D(q), Qij{q) the element oiQ{q) in Eq. (25). For steady state, Eq. (12), V-js(q) = 0. 



One can calculate that when Q = then jg = 0; if Q{q) ^ then generally js(q) 7^ 0, 
since dj [(5y(q)ps(q)] 7^ 0. Therefore the framework encompasses the cases without detailed 
balance. The term "detailed balance" used here is defined for abstract phase space, that is, 
the net current between any two states in the phase space is zero [31], identical to that for 
Markov process in mathematics. 



There are various ways to take the zero mass hmit (ehminate the fast degrees p), a 
derivation adopting the standard projection operator method [Tj is given in [TD], that is, to 
project the (q, p) space to the q space. We note that another commonly used approach of 
degree reduction, adiabatic ehmination [HlllISS], is based on the postulation that the velocity 
is a fast variable, which is consistent with the zero mass limit. Three other derivations 



for 1-dimensional systems are discussed in [TT]. This FPE (Eq. (25)) defines the A- type 



integration for SDK of the form of Eq. ([s]). Since generally Eq. ([I]) is equivalent to Eq. (|3 



Eq. (25) defines a stochastic integration for general Langevin equations. The steady state 



distribution in the phase space (if it exists) is the Boltzmann-Gibbs distribution Eq. (12). 
We should clarify that it is not to say that A-type is the only way to integrate Eq. (|3]), but 
a stochastic interpretation with physical significance. 

III. DISCUSSIONS 

A. Physical Meaning of Equation ([s]) 

We can consider a time-evolution system described by SDE as a particle motion inside 
the phase space. From a physical point of view, the motion of an object should be driven by 
some underlying forces. We consider this particle as a charged and massless one driven by 
the force: Ffotai = ""^q = (m = 0). The driven force can be separated generally in physics 
into three parts, dissipative, conservative and random: 



total " dissipative 



~r -T conservative ~r -T random U . v^ ' / 



Equation ^ has exactly this structure. In the usual study of 2 and 3 dimensional cases, 
Eq. ([3]) corresponds to a known fundamental equation in physics. We can generally use a fric- 
tional force as the dissipative component Fdissipative = — 'S'(q)q where S{q) is symmetric and 
positive semi-definite; a Lorentz force together with an energy induced force (for example, 
the electrostatic force) as the conservative component Fconservative = eq x B(q) + [— V0(q)]; 
the random component is described by a Gaussian white noise Frandom = ^{q)^(t) with zero 
mean that has the common origin with the frictional force, as formulated by the fluctuation 
dissipation theorem [36H38] : iV(q)iV^(q) = 2eS'(q), hence 

-5(q)q + eq X B(q) - V</.(q) + iV(q)e(t) = . (28) 
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The friction matrix S'(q) is positive semi-definite (guaranteed by the fluctuation-dissipation 
theorem), keeping the resistant property of the dissipative force (the allowed values are con- 
fined in the negative half space). It can be a non-diagonal matrix, describing an anisotropic 
frictional force. The potential function 0(q) plays the similar role to the Hamiltonian in 
dissipative systems, leading to the Boltzmann-Gibbs distribution on the final steady state (if 
it exists) of the stochastic process. For higher dimension, we generalize B(q) x q as A(q)q, 
where ^(q) is an antisymmetric matrix. There is a direct and simple physical realization 
based on this intuitive explanation of the SDE Eq. (|3|, for instance, practical experiments 
to verify the generalized Einstein relation are designed in Section 5.3 of [LS^. A recent 
experiment p^ has been implemented under a similar physical setting. 



B. Difference and Relation of the A-type Integration with Traditional Ones 

We have mentioned that the A-type integration reduces to an a-type one with a = 1 for 
one dimensional system, hence it is different from Ito or Stratonovich's. A more detailed 
demonstration on this point is provided in fll] . For general cases, when dimension is higher 
than one, things become more complicated, A-type integration is even not a-type, that 
means at different time intervals the choices of a are not the same, or the parameter a is 
position dependent a(q). A relation between A-type integration and Ito's is formulated in 
[16] by considering the correspondence of their FPEs, that is: the SDE 

q=(f + Af)(q) + iV(q)e(t), (29) 

where Af(q) = (A/i(q), . . . , A/„(q))^ using Ito integration describes the same process to 
that of the SDE 

q = f(q) + iV(q)e(t), (30) 

using A-type integration. Here A/j(q) = edj[Dij{q) + Qij{q)]. We note that there is an 
obvious distinction between Af (q) here and the additional (to Ito's) drift term for a-type 
(a choosing as a constant) integrations A/-- (q) = a [dkNij{q)] Nkj{q). 

An illustrative example has been provided in [16] to show a significant difference on 
the final steady state distributions between Ito and A-type integrations, the dynamics q = 

11 



f(q) + iV(q)eWis[39] 

X = 2x-x{x'' + y^) + ^2{x^ + y^)Ut) ,^^. 

y= 2y~y{x' + y^) + y/2{x^ + y^)^y{t) 

with q = {x, yY, e = 1, the diffusion matrix D{x, y) = {x^ + y'^) ■ I2, ^(t) = (^j,., C^yY has zero 
mean, and {C,{t)C,^{t')) = 6{t — t')l2- The final steady state distribution of this system with 
A-type integration can be calculated based on the framework as Psa{.x^ y) = ~zl7i ^-^P [~*^/^] 
with 

0(q) = -In [x' + y') + ^^ . (32) 



The distribution of Ito integration for the system Eq. (31) is identical to the A-type's dis- 
tribution for the system: 

' X = -x{x^ + y^) + ^2{x^ + y^)Ut) .33. 

y = -y (x' + y^) + ^2{x^ + y^Yyit) 

whose expression can be similarly calculated as psi{x,y) = ^^exp [—ip/e] with 

H^) = l{x' + y') . (34) 

Similarly, Ito integration for the system 

x = Ax-x{x^ + 2/2) + ^2{x^ + y^)Ut) .3^. 

y = Ay-y{x' + y") + ^ 2{x' ^ y'')iy{t) 

will reach the distribution psA^x^y). 

The two distributions psA and psi have obvious differences, for instance, Psa(0, 0) = 
but the origin (0, 0) is the most probable state for psi- These theoretical results have been 
verified by numerical experiments in [16] (shown in Fig. II]). Note that the potential 0(q) 



in Eq. (32) serves as a global Lyapunov function [17j for the deterministic dynamics (when 



e = 0) of the system Eq. (31). The Lie derivative [10] of 0(q) (V0(q) ■ f(q)) keeps non- 
positive. Therefore, there is an exact correspondence between the deterministic dynamics 
and the steady state distribution derived using A-type integration for SDEs, for example, 
the stable fixed points are locally most probable states, in accordance with one's intuition. 
This direct correspondence can not be kept by applying Ito nor Stratonovich integration: 
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(a)It6 Integration 



(b)A-type Integration 



FIG. 1. Sampling Results: The figure shows the distributions in the phase space for the long-term 



sampling of Eq. (31 ) using Ito and A-type integrations, the basal plane denotes the two dimensional 



phase space, the vertical axis indicates the emerging probability density of the system at a specific 
state. Fig. (a) is for Ito integration, it is consistent with the theoretical result Psi{x, y) oc exp [—ip/e] 



{il) is given in Eq. (34), with (0, 0) being the most probable state). Fig. (b) is for A-type integration, 
it is coherent with the deterministic dynamics obtained when e = 0, for example, being most 
probable at the circle x^ -|-y^ = 2 (are stable fixed points when e = 0), and having zero probability 
at (0, 0) (is an unstable fixed point when e = 0). 



After using Ito's, the unstable point (0,0) of Eq. (31) (e = 0) becomes stable in '?/'(q); 
moreover, stable points on x^ + y^ = 2 in Eq. (31) (e = 0) disappear in Eq. (34). The unique 



advantage of using A-type integration comparing to Ito or Stratonovich's is shown here. It 
enables a straightforward calculation of the transition probability from one fixed point to 
another after taking into account the noise infiuence for ODE models which is essential in 
many applications [3 [TJ El [121 El- For Ito or Stratonovich integrations, this is not direct 
and can even be impossible, since one can not recognize the original stable fixed points from 
the long time sampling distribution, like what shows in this example. 



C. 1-dimensional Example 

A comprehensive discussion on 1-dimensional systems has been conducted previously 
[llj . The corresponding FPEs and the steady state distributions for three different type of 
stochastic integrations (contained in [IT]) are listed in Table. |l| As has been shown in |¥T] . 
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there exist real processes choosing Stratonovich integration. A recent experiment [18] on 
a one-dimensional physical process suggested that A-type integration of the corresponding 
SDE is consistent with the experimental data. We show in the following that the zero mass 
limit is well established based on the physical setting of the experiment. 

For the colloidal particle studied in [TS], p = 1510 kg/m^, R = 655 x 10~^ m, then the 
volume V = A/SttR^ = 4/3 x tt x (655 x 10"^)^ m^ = 1.18 x 10"^^ m^, the mass of the 
particle is m = pV = 1510 x 1.18 x 10"^^ kg = 1.78 x 10"^^ kg. The symbol dt denotes the 
sampling time interval in [18], within dt < 10 ms, the authors claim that "the force acting 
on the particle can be treated as locally constant." In different experiments, they choose dt 
around the magnitude of 1 ms. 

We use the simplified formula for the friction coefficient from Mannella et a/.'s discussion 
of the same experiment [1], with r] = 8.5 mPa s and zq = 700 nm: 



-f{z) = 67rr]R^-^^^ . (36) 



For 2; — )■ oo, we have 700 = GirriR = 67r x 8.5 x 10 ^ x 655 x 10 ^ Ns/m = 1.1 x 10 '^Ns/m, 
then the time needed for the particle to return to equilibrium after acted by an external force 
F{z) is t{z) = [F{z)/-f{z)]/[F{z)/m] = m/-f{z), t^ = m/-f^ = 1.78 x lO^^Vl.l x 10"^ s = 
1.6 X 10"® s « 1 X 10"^ s = 1 ms ^ dt. In Fig. 2 of [18], for instance, at z = 200 nm, 
7(;z) = 1.1 X 10"^ Ns/m X (200 + 700)/200 = 5.0 x 10"^ Ns/m. Since -f{z) is always 
larger than 700, such that t{z) < too << dt. We conclude that during the sampling time 
interval the particle has been in equilibrium, hence the mass can be considered as zero in 
the experiment, realizing the zero mass limit. 

TABLE I. Fokker-Planck equations and steady state distributions with e = 1 
a Fokker-Plank equation Steady state distribution 



mP(<l^t)=[i,ijDiq) + l^D{q)^iq)]p{q,t) ^ |^ exp [-c/>(g)]} 



1 d_ 

2 dtf 



p{q,t)= lDkq)lDkq) + lD{q)(P{q) p{q,t) l -l exp [-0(g)] 



Zs 



D2(g) 



|p(g,t)= lD{q)i- + i-D{q)ct>{q) p{q,t) 1 exp [-0(g)] 



at' 



' dq dq 
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D. Previous Attempts 

Some related previous attempts have been discussed in our former works [H |9l [13] . In 
this subsection we provide a brief summary. No effort is made to have a complete list, while 
we do hope to have picked some of the major ones. 

To our knowledge, the one dimensional case of the present framework, a = 1 type inte- 
gration, has been mathematically discussed by Wong and Zakai in the 1960s [12j (see also 
P [13]). Their work provided a method for a G [0, 1]. 

From a physics perspective Grabert et al. |13] provide a method for nonlinear irreversible 
processes, and was discussed by Graham in a long review which we had commented [6] . Their 
method has a restriction based on the singularity of diffusion matrix -D(q) (see Eqs. (2.14- 
16) in |l3]) and handles the phase variables separately according to a reversibility defined 
by such singularity. By contrast, our framework can be used directly for a given -D(q) in 
its general form without a separated treatment, a global potential is then obtained in the 
whole phase space. 

Theoretical considerations for the situations with detailed balance have been proposed by 
Klimontovich [44j (see also [13j), while our approach applies to situations of both with and 
without detailed balance. On the weak noise limit, previous works reach a same Hamilton- 



Jacobi equation to our Eq. (11), the nondifferentiable potentials reported in [l5l[l6] are then 
included within our framework; moreover, our potential function is valid for an arbitrary 
noise strength. A major difference between our construction of the potential and those in the 
literature such as Graham-Haken construction |17] is that no assumption on the stationary 
distribution function is needed in our approach in the limit t — > oo. In particular, the 
potential in our framework can be time dependent (see also [H [20]). Other results beyond 
Ito and Stratonovich have been proposed in [18H50] . 

IV. CONCLUSION 

Beginning with a 2n dimensional stochastic differential equation (SDK) system we de- 
rived a generalized Klein-Kramers equation. It is a Fokker-Planck equation (FPE) that can 
be obtained from the SDK system regardless the type of stochastic calculus used (Ito and 
Stratonovich integrations of this SDK are described by the same FPE). This provides a nat- 
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ural starting point to demonstrate the internal consistency of our approach. After taking the 
zero mass hmit, the 2n dimensional SDE reduces to an n dimensional structured form that 
is equivalent to the n dimensional Langevin equation with multiplicative noise based on our 
assumption. The corresponding generalized Klein-Kramers equation turns to a new FPE. 
Such limiting process defines a new type of stochastic integration (A-type) for Langevin 
equations different from traditional ones. During the demonstration of consistency, the Ito 
vs. Stratonovich dilemma is not encountered. In addition, the new FPE is generally not 
reachable by the a-type integration in higher dimensions. The framework and the A-type 
stochastic integration are consequences under a physical view of time evolution dynamics 
in the phase space. An attractive advantage of A-type integration is the natural correspon- 
dence between stochastic and deterministic dynamics. For example, fixed points are not 
changed, which is not possessed by Ito or Stratonovich's. As an illustration, experimental 
data demonstrates that a one-dimensional physical process realizes the zero mass limit and 
chooses A-type integration. 
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